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Abstract 

We reexamine constraints on the spatial size of closed toroidal models with 
cold dark matter and the cosmological constant from cosmic microwave back- 
ground. We carry out Bayesian analyses using the Cosmic Background Ex- 
plorer (COBE) data properly taking into account the statistically anisotropic 
correlation, i.e., off-diagonal elements in the covariance. We find that the 
COBE constraint becomes more stringent in comparison with that using only 
the angular power spectrum, if the likelihood is marginalized over the orien- 
tation of the observer. For some limited choices of orientations, the fit to the 
COBE data is considerably better than that of the infinite counterpart. The 
best-fit matter normalization is increased because of large-angle suppression 
in the power and the global anisotropy of the temperature fluctuations. We 
also study several deformed closed toroidal models in which the fundamental 
cell is described by a rectangular box. In contrast to the cubic models, the 
large-angle power can be enhanced in comparison with the infinite counter- 
parts if the cell is sufficiently squashed in a certain direction. It turns out that 
constraints on some slightly deformed models are less stringent. We comment 
on how these results affect our understanding of the global topology of our 
universe. 
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I. INTRODUCTION 



Is the universe finite or infinite? If we assume the "strong Copernican principle" (SCP) 
that the spatial hypersurface is globally and locally homogeneous and isotropic, then all 
three-spaces are classified by the sign of the curvature; spaces with a positive constant 
curvature are closed and finite in volume while spaces with a negative constant curvature 
or vanishing curvature are open and infinite. Hence, one might expect that the sign of the 
spatial curvature completely determines the size of our universe. 

However, there is no particular theoretical reason for assuming the global homogeneity 
and isotropy of the spatial hypersurface, since the Einstein equation itself cannot specify the 
global topology of the space-time. In fact, there are a variety of cosmological models that 
satisfy the "weak Copernican principle" for which the spatial section is locally homogeneous 
and isotropic but not necessarily globally homogeneous and isotropic. These multiply con- 
nected (MC) models (for a comprehensive review see [1-3]) are topologically distinct but 
the metric is exactly the same as that of the simply connected counterparts with constant 
curvature. Therefore, these MC models are in harmony with the observational facts that 
support the Friedmann-Robertson- Walker (FRW) models, namely, the cosmic expansion, 
the relative amounts of primordial light elements, and the cosmic microwave background. 
In order to answer whether or not the spatial section of our universe is finite, we need to 
determine not only the curvature, but also the global spatial topology of the spatial hyper- 
surface, since there are a number of closed compact spaces with a negative or a vanishing 
curvature. 

If the spatial section is sufficiently small, then all the fluctuations beyond the comoving 
size of the space are suppressed. It has been claimed by several authors that the "small 
universes" [4] in which the topological identification scale is equivalent to or smaller than 
the horizon scale have been observationally ruled out because large-angle power of the cosmic 
microwave background (CMB) temperature fluctuations is significantly suppressed [5-8] and 
does not fit to the Cosmic Microwave Background Explorer (COBE) Differential Microwave 
Radiometer (DMR) data. 

However, for low matter density models Q < 1, constraints on the size of the spatial 
section are less stringent, since the bulk of large-angle temperature fluctuations can be 
produced in the curvature or A dominant era, leading to less stringent suppression in the 
large-angle power [9-12]. 

On the other hand, it has been pointed out that constraining closed MC models using 
only the angular power spectrum is not sufficient, since it does not contain information 
about the anisotropic component of statistically averaged fluctuations [13]. In fact, temper- 
ature fluctuations in the sky for the MC models (except for the projective three-space RP 3 ) 
are described by an anisotropic Gaussian random field for a particular choice of position 
and orientation of the observer, provided that each Fourier mode is independent and obeys 
Gaussian statistics as predicted by the standard inflationary scenario. Therefore, in order 
to constrain the MC models, one must compare all anisotropic and inhomogeneous fluctua- 
tion patterns for every possible choice of orientation and position of the observer with the 
observed data, which needs time-consuming Monte Carlo simulations. Based on the study 
of several closed hyperbolic (CH) models, [13] claimed that all "small universes" are ruled 
out because the anisotropic correlation patterns are at variance with the COBE DMR data. 
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Nevertheless, for CH models, a subsequent analysis confirmed that the fit to the data 
depends not only on the orientation but also on the position of the observer. In some limited 
places, the likelihood is considerably higher than that of the simply connected counterpart 
[14] . Even if the place at which the pixel-pixel two-point correlation agrees with the data is 
limited, the model cannot be ruled out if the fit to the data is sufficiently good at certain 
places. 

On the observational side, there has been mounting evidence that the spatial section of 
our universe is almost flat [15-18]. If the background space is exactly flat, the topology of 
the orientable closed space is limited to six kinds [1]. The simplest space is the toroidal 
space T 3 . Although closed flat MC models need a "fine-tuning" that the present horizon 
scale is comparable to the topological identification scale Lt, 1 it is of crucial importance to 
give a lower bound on the size of the spatial section of these models, since we have no prior 
knowledge of the spatial topology of the universe. 

In this paper, we improve the previous COBE bounds on the size of the spatial hypersur- 
face of closed toroidal flat cold dark matter models with a consmological constant (ACDM) 
[19] by carrying out Bayesian analysis properly taking account of statistically anisotropic 
correlation that cannot be described by the angular power spectrum C/. The temperature 
fluctuations are calculated by solving the linear evolution of the photon-baryon fluid from the 
COBE scale to the acoustic oscillation scale. We estimate the best-fit normalization of the 
power of the matter density fluctuations at the 8h _1 Mpc scale a 8 in comparison with those of 
simply connected ACDM models. Throughout this paper we assume that Q m = 0.3, Q\ = 0.7. 
In Sec. II, we give a brief account of the temperature anisotropy on large to intermediate 
scales and we study how the size and the shape of the fundamental cell affect the large-angle 
anisotropy. In Sec. Ill, we carry out full Bayesian analyses using the COBE DMR data based 
on pixel-pixel correlation and obtain the likelihoods and the best-fit normalizations a$. In 
Sec. IV, we discuss the viability of other closed flat models, almost-flat closed spherical and 
almost-flat closed hyperbolic models. In Sec.V, we draw our conclusions and discuss the 
possibility of detecting the signature of the finiteness of the spatial geometry of the universe 
in the future. 



II. TEMPERATURE ANISOTROPY 



In the Newtonian gauge, the /th multipole of the primary temperature anisotropy in 
the FRW models projected onto the sky at the conformal time 77 is written in terms of the 
monopole ©0, and the dipole ©1 of temperature fluctuation and the Newtonian potential ^ 
at the last scattering and the time evolution of ^ and the Newtonian curvature $ 



= [©o + k)mx(v) - xfa.))) + ©i(^, - xfaO)) 



In this paper, we define L t as the comoving length of a closed geodesic which connects a point x 
and another point gx in the universal covering space where g is a generator of the discrete isometry 
group T that defines a Dirichlet fundamental domain (i.e., a fundamental cell). 
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+ £(^-^>(*(x(»7)-x(»/)))d»/, (1) 

where k denotes the wave number, 77* is the last scattering conformal time, the radial distance 
is xiv) = V an d Ji represents the /th-order spherical Bessel function. The first term in the 
right-hand side in equation (1) represents the "ordinary" Sachs- Wolfe (OSW) effect caused 
by the fluctuations in the gravitational potential at the last scattering. The second term 
describes the Doppler shift owing to the bulk velocity of the photon-baryon fluid. Note that 
this term is negligible at superhorizon scales kq < 1 since radiation pressure cannot play any 
role and cosmic expansion damps away the initial velocity. The third term represents the 
integrated Sachs- Wolfe (ISW) effect caused by the decay of the gravitational potential. 

In order to constrain the FRW models that yield isotropic Gaussian fluctuations, one 
needs only two-point correlations of temperature fluctuations AT/T in the sky, which can 
be characterized by the angular power spectrum Ci = (|a^ 12 v 



tin 



^' c , ' r* t M, (2) 



4tt 2tt 2 Jo k 2/ + 1 

where a\ m represents the expansion coefficient with respect to the spherical harmonic Yi m , 
() denotes an ensemble average over the initial fluctuations, and rjo is the present conformal 
time. Note that C\ is invariant under any 5*0(3) transformations. 

Now, let us consider a closed toroidal model T 3 in which the spatial hypersurface is 
described by a cube with side L in the Euclidean three-space whose opposite faces are 
identified by translations. Then the wave numbers of the square-integrable eigenmodes of 
the Laplace-Beltrami operator are restricted to the discrete values ki = 2nni/L, {i = 1, 2, 3) 
where n^s run over all integers. Equation (2) now reads 

2i + i„ 1 ^ \e,(m,k)\ 2 ,,, 

-s-^E-ann- p) 

Note that all the modes whose wavelength is longer than L are not where ai m represents 
the expansioallowed ("mode cutoff"). From now on, we restrict ourselves to the adiabatic 
scale-invariant Harrison-Zel'dovich spectrum n — 1 [i.e., V${k) = /c 3 (|$(0, k)\ 2 ) = const.] as 
predicted by the standard inflationary scenario. Note that ATJ = [l(l+l)Ci/ (2tt)} 1 ^ 2 becomes 
almost flat for the Einstein-de Sitter universe. For low-density FRW models, the large-angle 
power is boosted owing to the ISW effect, since A dominant epochs come earlier. However, 
for T 3 models, the OSW contribution to the large-angle power is suppressed because all 
modes whose wavelength is longer than the cutoff scale L are not allowed. The angular 
scale below which the power is suppressed is approximately given by l cut = 27ir] /L — 1 [19]. 
Interestingly, for low-density T 3 models with moderate size of the cell, the excess power 
owing to the ISW contribution is almost canceled out by the suppression owing to the mode 
cutoff. As shown in figure (1), the suppression in the large-angle power is not prominent 
even for a model with e = L/2rj = 0.2. 
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FIG. 1. Plots of the angular power spectrum AT] = [1(1 + 1)Q/(2tt)] 1/2 for cubic toroidal 
models with various linear sizes parametrized by the ratio of the length of the side to the diameter 
of the last scattering surface in the comoving coordinates :e = L/2rjQ in comparison with that for the 
simply connected counterpart (solid curve). AT; is normalized by ATio for the simply connected 
models. Cosmological parameters are $7a = 0.7, Qcdm = 0.26, Of, = 0.04, and h = 0.7. 
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FIG. 2. Plots of the angular power spectrum AT} = [(/ + l)Ci/ {2-k)] 1 / 2 for deformed toroidal 
models whose fundamental cell is described by a rectangular box with sides L\=L2, L3. The volume 
is fixed to V = (0.4 x 2t]q) 3 . The degree of deformation is parametrized by s = (L\ = L^jLj,. 
ATi is normalized by AT10 for the simply connected counterpart. In the small- angle limit I > 1, 
the angular power spectrum converges to that of the simply connected counterpart. Cosmological 
parameters are the same as in figure 1. 



Next, we consider models in which the fundamental cell is described by a rectangular 
box D with length Lj of sides. We call these models the "pancake' type or the "cigar" type 
if D has a relatively long side in one direction, or long sides in two orthogonal directions, 
respectively. As shown in figure 2, the large-angle power is enhanced for models with a 
sufficiently squashed cell. It seems that the enhancement is rather surprising since Fourier 
modes whose wavelength is larger than the size of the cell are not allowed. The apparent 
enhancement is related to the number density of the Fourier modes. In order to see this, 
let us consider the number function defined as the number of eigenmodes of the Laplace- 
Beltrami operator whose eigenvalue is less than E = k 2 , which is described by Weyl's formula 
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in the small-scale limit N(E) = Vk 3 /27r 2 , where V denotes the volume of the three-space. 
For globally "almost" isotropic spaces like cubic toroidal spaces, the Weyl's formula gives a 
good approximation of N(k) in the large-scale limit as well. However, for globally "very" 
anisotropic spaces 2 , the deviation from the Weyl formula is prominent in the large scale 
limit. For instance, the number function is approximately given by N oc k for the "cigar" 
type and iV oc k 2 for the "pancake" type [21]. The discrete wave number is given by 



k = 2tt, 



2 2 2 

Hi + Hi j_ 3 



2 ' 

3 



(4) 



where each n« runs over all integers. Hence, for L 2 ,Ls <C L±, the modes on large scales 
k <C 1/L 2 ,1/L 3 are given by k = 2im\/Li, leading to N(k) ~ Lik/n. Similarly, if L 3 <C 



Li ~ L 2 , the modes on large scales k <C L 3 are given by k = 2tt^^ + ^|, then N(k) ~ 

LiL 2 k 2 /4tt. Because the power is normalized by the total volume, the contribution of modes 
on large scales are relatively boosted, resulting in an enhancement in the large-angle power. 
Interestingly, for slightly squashed models, the reduction in the large-angle power owing to 
the mode cutoff is almost canceled out by the enhancement owing to the squeezing of the 
space. 

Although we have assumed that the primordial power spectrum is scale-invariant even 
in the case of deformed models, the non-trivial boundary condition on the mode functions 
of inflaton fields may cause a deviation from the scale-invariant spectrum. However, such 
a deviation is expected to be small, since the background space time has a high degree of 
symmetry locally at the inflationary epoch that is relevant to the n — 1 spectrum. 

So far, we have studied the effect of the nontrivial global topology on the large-angle 
power. However, the angular power spectrum is not sufficient for constraining the models 
since the temperature fluctuations in the MC models are not statistically isotropic in general; 
the effect of off-diagonal elements (aj m a*/ m /) for I ^ V or m ^ m' which are not SO (3) 
invariants cannot be negligible. Therefore, it is necessary to perform Bayesian analyses 
taking account of these off-diagonal elements for each possible orientation of the observer, 
since the toroidal models are globally anisotropic. 

III. BAYESIAN ANALYSIS 

Any Gaussian temperature fluctuations in the sky can be characterized by a covariance 
matrix whose elements at pixel i and pixel j are, 

Mij = {T % T 3 ) = (a lm 4 m ,)W l W l ,Y l Un l )Y l , m ,(h :i ) + (AW,) (5) 

ll'mm' 

where ai m is an expansion coefficient with respect to a spherical harmonic Yj m , and () 
denotes an ensemble average taken over all initial fluctuations for a fixed orientation of the 



2 The degree of global anisotropy in the geometry can be measured by the ratio of the maximum 
geodesic shortest distance between two arbitrary points to the characteristic length scale V 1 / 3 . 
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observer. Tj represents the temperature at pixel i, Wf is the experimental window function 
that includes the effect of beam smoothing and finite pixel size, rij denotes the unit vector 
toward the center of pixel i, and (iVjJVj) represents the noise covariance between pixel i and 
pixel j. 

If we assume a uniform prior distribution for the cosmological parameters, from Bayes' 
theorem, the probability distribution function of a set of C\™ = (ai m a* Um ,) for a random 
Gaussian field (not necessarily isotropic) is given by 

A(C '™ ' |f » K def4 (C£ -) 6XP fi fT ' M_1(CL "" ) ' f ) ■ <6) 

In the following analysis, we use the inverse-noise-variance-weighted average map of the 
53A, 53B, 90A, and 90B COBE DMR channels. To remove the emission from the galactic 
plane, we use the extended galactic cut (in galactic coordinates) [22]. After the galactic cut, 
the best-fit monopole and dipole are removed using the least-squares method. To achieve 
efficient analysis in computation, we compress the data at "resolution 6" (2.6°) 2 (I ~ 72) 
pixels into one at "resolution 4" (10. 4°) 2 (Z ~ 18) pixels for which there are 384 pixels in the 
celestial sphere and 248 pixels surviving the extended galactic cut. The window function is 
given by W\ = GiF t where F t are the Legendre coefficients for the DMR beam pattern [23] 
and Gi are the Legendre coefficients for a circular top-hat function with area equal to the 
pixel area, which account for the pixel smoothing effect [24] . We also assume that the noise 
in the pixels is uncorrelated from pixel to pixel, which is known to be a good approximation 
[25]. 




0.2 0.3 0.4 0.5 0.6 



FIG. 3. COBE constraints on cubic models. The solid curve represents the likelihoods A 
marginalized over 24000 orientations of the observer and the normalization for various relative 
linear size e = L/2t]q. The dashed line shows the maximum values. The dotted line represents the 
likelihoods computed by neglecting the anisotropic components. All the likelihoods are normalized 
by that of the simply connected counterpart A(oo). Cosmological parameters are the same as in 
figure 1. 
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Because the toroidal T 3 model is globally anisotropic, the two-point pixel-pixel correla- 
tions depend on the orientation of the coordinate axes. Therefore, we need to compute the 
likelihood function for each orientation a of the observer. Assuming a constant distribution 
for a and the normalization a, the likelihood is given by A = / A(a, a)dada, where da 
denotes the volume element of a Lie group SO (3) with Haar measure. 

As shown in figure 3, the angular power spectrum alone does not give a stringent con- 
straint e > 0.3 on the size of the cell. However, if we use a likelihood function taking account 
of anisotropic two-point correlations marginalized over the orientation of the observer, we 
obtain a rather stringent constraint. 

Because there is no upper limit for e, we choose a variable N = n/6e 3 = (volume of the 
observable region at present) /(volume of space at present) as a physical quantity that is 
assumed to be constantly distributed. We find that the logarithm of relative likelihoods 
log 10 [A(A)/A(A = 0)] for cubic models (0.2 < e < 0.65) can be well fitted by a polynomial 
function —bN a where a = 0.46, 6 = 0.78 which are obtained by minimizing \ 2 - By marginal- 
izing over N, we find a constraint N < 2.4 (e > 0.6) at 68% CL, and N < 8.2 (e > 0.4) at 
95% CL corresponding to log 10 A/A(oo) ~ —0.9, log 10 A/A(oo) ~ —2.0, respectively. The 
constraint on the size of the spatial section is more stringent compared with the previous 
analysis [19] in which fluctuations on smaller angular scales l> 10 have not been taken into 
account. 

It should be emphasized that the likelihoods depend sensitively on the orientation of the 
observer, especially for models with small sides e < 1. The maximum value of the likelihood 
far exceeds the mean by a factor A max /A = io 2 - 5 ~ 3 - 9 for 0.2 < e < 0.7 which suggests that 
almost all the orientations are ruled out. For instance, A max /A ~ 10 3 ' 9 implies that the 
likelihood attains high values only when the axes are oriented with a precision of ~ 5° which 
is equal to one half the size of one pixel. 3 

To see the orientations that give the enormous likelihood, we have plotted the orientations 
in the skymap that give a large likelihood. Except for e = 0.2, in which a weak correlation 
to the galactic pole is found, no particular preferred orientation is found (figure 4). 



3 For 24000 realizatons of orientation, the distance between a pair of the nearest two points on 
SO(3) is approximately a cubic root of [volume of SO(3)]/24000/(number of symmetries of T 3 ) = 
(87r 2 /24000/24) 1//3 x 180/7T ~ 3°. Therefore, each peak (~5°) is represented by just several points. 
We have checked the convergency of the calculations by increasing the number of points near the 
peaks to several hundreds, which resulted in just 10 percent difference in the relative likelihoods. 
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FIG. 4. Best-fit orientations in the skymap for five cubic models. Each orientation is described 
by a set of six points that correspond to directions to the centers of the face of the cube that are 
preserved by the symmetries of T 3 . A total of 40 orientations that give more than 80 percent 
contribution to the averaged likelihoods are plotted for each model. All the orientations for the 
five models are shown together in the last figure. The Aitoff projection is used. The center of the 
oval represents the galactic center. 




FIG. 5. The co ntribution of off-dia gonal elements for cubic models. The plotted lines repre- 
sent | (ai m a*, m , )\/^J {ai m a* Vm , ) (ai m a*, m , ) for 2 < I < 15, in ascending order of 1(1 + 1) + m + 1 and 
I' (I' + 1) + m' + 1 for two types of cubic models. 
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In the previous literature, the effect of the global anisotropy has been considered to be 
small for cubic cases, because the contribution of off-diagonal elements is quite small. This 
argument is valid as long as one considers the averaged values over all the possible set of 
parameters (/,/', m, ml). Let us consider a set of orthogonal coordinates (x,y,z) in which 
all the faces of the cube are perpendicular to one of the coordinate axes. Because of the 
symmetries of T 3 , (ai )Tn aJ, m ,) vanishes if / + /' is even and (m + m')/2 is even. In fact, the 
contribution of off-diagonal elements relative to the diagonal elements are small as 0.042 to 
0.015 for e =0.2 to 0.6 for an ensemble of all values of (l,l',m,m'). However, as shown in 
figure 5, each nonzero component of off-diagonal elements is not so small (0.33 to 0.12 for 
e = 0.2 to 0.6). Although any rotation of the coordinates gives a different structure in the 
off-diagonal components, the spiky structure cannot be completely vanished by any rotation. 

One might argue that computing a likelihood marginalized over every orientation is 
meaningless and the constraint should be given by the maximum value of the likelihood 
instead of the marginalized value. Then the constraint becomes less stringent e > 0.25. 
However, if we assume that the universe does not prefer a specific orientation, the cases 
of the "failed orientations" should be taken into account. By deforming the shape of the 
fundamental cell, or by changing the amplitude of the fluctuations, there might be a chance 
that the fluctuation fits the data almost perfectly, since the data themselves are anisotropic. 
In order to check this, a similar analysis has been carried out for deformed toroidal models 
whose fundamental cell is described by a rectangular box with sides L 1 = L 2 and L 3 . As 
shown in figure 6, the "pancake" type models {L\ = L 2 )/L 3 > 1 give much better fits to 
the data. For instance, for s = 10 1 / 2 , the likelihood is improved by 10 times in comparison 
with the cubic model with the same volume V = (0.4 x 2r]o) 3 . The corresponding number of 
copies of the fundamental cell inside the last scattering surface in the comoving coordinates 
is N = 8 and the minimum linear size of the cell is equal to L 3 /2r] = 0.19. 
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FIG. 6. COBE constraints on deformed models whose fundamental cell is described by a rect- 
angular box with sides L\ = L2 and L3. The volume is fixed to V = (0.4 x 2r]o) 3 . The solid 
curve represents the likelihoods A marginalized over 24000 orientations of the observer and the 
normalization for various models. The dashed line shows the maximum values. The dotted line 
represents the likelihoods computed by neglecting the anisotropic components. All the likelihoods 
are normalized by that of the simply connected counterpart A(oo). Cosmological parameters are 
the same as in figure 1. 



In the previous literature, it has also been argued that an increase in the best-fit normal- 
ization constant is related to the large-angle suppression owing to the finite size of the spatial 
section [12]. From our analyses, it is found that inclusion of off-diagonal term can also give a 
comparable increase in the best-fit normalization constant as well. For instance, the best-fit 
normalization <j 8 is increased by ~ 10 percent for e = 0.6 (figure 7). The effect is relevant 
to the globally anisotropic geometry of the background space (see the appendix) which has 
been found for closed hyperbolic models (although the effect of the global inhomogeneity is 
neglected) [13]. 

If the best-fit normalization is too large, the choice of orientation is less probable (figure 
8) (see the appendix). We have checked that the anticorrelation between the likelihood and 
the best-fit normalization becomes weaker and the plotted points get more aligned as the 
comoving size of the spatial section increases. 
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FIG. 7. Plots of the ratio of the best-fit normalization of the cubic models to that of the simply 
connected counterpart (solid curve) and those obtained by an analysis in which the covariance 
matrix is computed by neglecting the off-diagonal terms (dotted curve). 
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FIG. 8. Plots of the best-fit normalization versus the likelihood for various orientations of the 
observer. The size of a cubic cell is e = 0.4. Each quantity is normalized by the value for the 
simply connected counterpart. 
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IV. OTHER TOPOLOGIES 



In addition to the toroidal models T 3 we have considered, there are five other 
flat orientable topologies, which can be realized as compact quotients of T 3 , namely, 
T 3 /Z 2 ,T 3 /Z 4 ,T 3 /Z 2 x Z 2 ,T 3 /Z 3 , and T 3 /Z 6 where Z m is a cyclic group of order m. The 
eigenmodes of the Laplace-Beltrami A operator can be written in terms of a finite sum of 
eigenmodes (= plane waves) of T 3 , 

«k( x ) K H expfik • (&x)], (7) 

where T denotes the discrete isometry group of T 3 and x is the Euclidean coordinates. Let 
us consider the case of T 3 /Z 2 , which can be realized by identifying the two pairs of opposite 
faces of a cube with sides L by translations but for one pair of faces with a translation with 
7r rotation. Because the finite-sheeted covering space is described by a rectangular box with 
sides L, L, 2L, one might expect that the space can support a mode with wavelength 2L. 
However, this is not true [26]. In fact, equation (7) yields the explicit form of the eigenmodes 

■Uk(x) oc exp[ik • x] + exp[—i(k 1 x 1 + k 2 x 2 + ik 3 (x 3 + L)], (8) 

where k\ = 2nn 1 /L,k 2 = 27rn 2 /L, k 3 = nn 3 /L. A mode with n\ = n 2 = 0,n 3 = 1 is 
equal to zero. Thus the allowed minimum wave number is still k = 2-k/L. Similarly, one 
can show that the maximum wavelength of the eigenmodes is comparable to the topological 
identification scale L t for other four cases [26]. Therefore, it is reasonable to conclude that 
the constraints on these closed flat models do not grossly differ from those on T 3 models. 

Although we have assumed that the background space is exactly flat so far, it should be 
emphasized that a slight deviation from the flat space means a significant difference in the 
global topology. For instance, if the background space is exactly flat, there is no particular 
scale that sets the limit on the comoving size of the topological identification scale L t . In 
other words, one can consider arbitrarily small or large spaces. However, for positively 
curved (spherical) or negatively curved (hyperbolic) spaces, the curvature radius R c sets a 
certain limit on the size of the space. 

For instance, it has been proved that closed hyperbolic spaces with volume larger than 
0.16668 • • ■ R z c cannot exist [27]. The smallest known CH space has volume V = 0.94R 3 [28]. 
Remembering that the volume of a ball with radius R (proper length) in hyperbolic space is 
given by V = 7r_R 3 [sinh(2_R/i? c ) — 2R/R C ], the topological identification scale of CH spaces 
with V ~ i? 3 is estimated as L t ~ R c . In this case, n~ 1/2 = Rc/H- 1 ~ Lt/H- 1 which gives 
L t ~ H- 1 for tt K = 0.1 and L t ~ lOH' 1 for tt K = 0.01. Thus, the effect of the global 
topology is almost negligible if fix ^ 0.1. However, if the spatial geometry is squashed 
in a certain direction while keeping the volume finite, the imprint of the global topology 
could be still prominent as in the case of flat topology 4 although the enhancement in the 
large-angle power may put a limit on the possible shape of the fundamental cell. Moreover, 



4 The fluctuations in "squashed" CH spaces resemble those in noncompact hyperbolic spaces with 
cusps [21]. 
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if we allow CH orbifold models with a set of fixed points, then the volume can be as small 
as V = 0.0391-R^ [14]. Then the effect of the nontrivial global topology is conspicuous even 
if the space is as nearly flat as Qk = 0.05 [14,29]. 

In the case of spherical topology, the volume is written in terms of the order |T| of the 
discrete isometry group T of SO(4), V = 27r 2 _R^/|r| [1]. Thus, the largest space is the 
three-sphere S 3 itself. On the other hand, there is no lower bound for the volume since 
one can consider a group with arbitrarily large number of elements of T. For instance, one 
can consider a cyclic group Z m with arbitrarily large m. Then V = 2n 2 Rl/\Z m \ can be 
arbitrarily small and the obtained quotient space is squashed in one direction. As in CH 
models, the imprint of the global topology can be prominent for squashed spaces, even if 
the space is as nearly flat as Vl K = 0.05 [30]. 

V. CONCLUSION 

In this paper, constraints on closed toroidal models from the COBE data have been 
reconsidered. By performing Bayesian analyses taking into account anisotropic correlation 
(i.e., off-diagonal elements in the covariance matrix) that has been often igonored in the 
previous literature, we obtained a constraint on the linear size of a cubic cell L > 0.6 x 2rj 
for Qa = 0.7, Q m = 0.3 based on the likelihoods that are marginalized over the orientation 
of the observer. The maximum allowed number N of copies of the cell inside the last 
scattering surface is ~ 2.4 (68% CL). We find that the constraint becomes more stringent in 
comparison with that using only the angular power spectrum. The best-fit normalization 
<t 8 can be increased by ~ 10 percent owing to the suppression in the large-angle power and 
the globally anisotropic geometry. We find that the likelihood depends sensitively on the 
orientation of the observer. Even if the likelihood that is marginalized over the orientation 
is small, the fit to the data for some limited set of orientations is far better than that of the 
infinite counterpart [14]. If one takes the maximum likelihood value instead, the constraint 
becomes conspicuously less stringent. Moreover, for some slightly deformed toroidal models, 
we find that the large-angle power can be flat and the constraint also becomes less stringent 
N ~ 8 (see also [31]). 

In order to detect the periodic structure owing to the nontrivial topology, the topological 
identification scale (more precisely twice the injectivity radius at the observing point) must 
be smaller than the diameter of the last scattering surface in comoving coordinates: L t /2r] < 
1. Therefore, the present constraints do not completely deny the possibility of detection of 
periodic structure by future observations [32]. 

Even in the case of L t /2r]o ~ 1, the signature of nontrivial topology might be observ- 
able as a nontrivial correlation structure in the temperature correlations [33] because the 
fluctuations in the MC models are non-Gaussian for isotropic and homogeneous observers. 
Inhomogeneous and anisotropic Gaussian fluctuations for a particular choice of position 
and orientation are regarded as non-Gaussian fluctuations for a homogeneous and isotropic 
ensemble of observers [34,35]. If non-Gaussian fluctuations with vanishing skewness but non- 
vanishing kurtosis are found only on large angular scales in the CMB, then it will be a strong 
sign of the nontrivial topology, or equivalently the breaking of the SCP [2,36-38]. Even if 
we failed to detect the identical patterns or objects in the sky, the imprint of "finiteness" 
could be still observable by measuring such statistical quantities. As we have discussed, the 
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MC spaces can be squeezed irrespective of the type of topology. Even if the volume of the 
spatial section is large, we might be able to detect globally anisotropic structure if the spatial 
section is sufficiently squeezed in a certain direction, The ongoing satellite missions such as 
MAP and Planck will provide us an answer to the question "how large is our universe?" . 
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APPENDIX A: ANISOTROPIC GAUSSIAN FLUCTUATIONS 



For homogeneous and isotropic Gaussian fluctuations, the variance of the expansion 
coefficient a\ m is constant in m for a given / regardless of the position and orientation of 
the observer. In other words, there is no "m-structure" or "direction-direction" correlation 
in the expansion coefficients. On the other hand, for anisotropic Gaussian fluctuations, the 
variance of a/ m depends on m [34]. Because a; m 's extracted from the data actually depend 
on m owing to the cosmic variance, for some very limited choices of the orientation, the 
"m-structure" in the anisotropic fluctuations may agree with the apparent "m-structure' in 
the data. 

To confirm this, we compare the variances of expansion coefficients b\ m of the best-fit case 
to those extracted from the data, which can be obtained by multiplying the temperature 
fluctuation T; L and the standard deviation in the noise iV(fij) by spherical harmonics Yj m , and 
carrying out an integration over the pixels surviving the "extended" galactic cut (denoted 
as SK) 

b lm = J2a Vm/ w v [ Y l , m ,(n)Y l i(n)dn+ f N(n)Y^(n)dn. (ai) 

I'm' JSK JSK 

As shown in figure 7, an excellent agreement at I — 6 — 7 is observed. Because the theoretical 
prediction of the quadrapole is somewhat higher than the observed value, the large-angle 
suppression owing to the finite size of the spatial section is not relevant. We also checked 
other cases that give a large value in the likelihood compared to that of the simply connected 
counterpart and again such an excellent agreement in \bi m \ at certain range of angular scales 
is observed for each case. Thus, for a limited choice of orientation, globally anisotropic 
Gaussian fluctuations can give a better fit compared with globally isotropic and homogeneous 
Gaussian fluctuations. 




FIG. 9. Plots of \bi m \ = yj (bi m b* m ) in ascending order of 1(1 + 1) +m + 4 for the T 3 cubic model 
(solid curve) with e = 0.4 (solid curve) which best fits bi m extracted from the COBE DMR data 
(dashed curve). 
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Next, we consider the effect of off-diagonal elements in the two-point correlation 
(a,im a i' m 1 ) i I ^ I' oi m ^ m! , which represent scale-scale correlations. 

For illustrative purpose, first, we consider a toy model in which the fluctuations obey a 
bivariate Gaussian distribution function, 



p(yi,y 2 ,cr, r) 



1 



cxp 



'vl + vt- ^y x y 2 



2(1 - r 2 



(A2) 



where r = cov(y 1 ,y 2 )/a 2 is the correlation coefficient which represents the off-diagonal 
components. By an orthogonal transformation, y\ = (xi + x 2 )/y/2 : y 2 = {—X\ + x 2 )/\/2, we 
obtain a distribution function written in terms of a diagonalized covariance matrix, 



p(x 1 ,x 2 ,a,r) 



27ra 2 V / l^ 2 ~ 



exp 



(l-r)x 2 + (l + r) 



~, 2 
J 2 



2(l-r 2 )a 2 



(A3) 



For a given set of X\ and x 2 , the best-fit normalization cr max is given by 



(1 - t)x\ + (1 + t)x\ 
2(1 -r 2 ) ; 



(A4) 



for which p(a) takes a maximum value. Suppose a distribution p with r = 0. A slight 
increase in r causes a squeezing of the distribution p in the direction x\. Then, in the 
region \x 2 \ > \x±\, o" max increases while in the region \x 2 \ < \x±\ o" max decreases. As r 
increases, the distribution is much squeezed and the region in which <r max increases widens 
as | #2 1 > \J(1 — t)/(1 + t). Therefore, in comparison with the case in which r = 0, for a 
wide range of parameter region (xi, x 2 ), <r max is increased if r is large enough. On the other 
hand, the maximum value of the probability decreases, since it is inversely proportional to 

Although this toy model is over simplified, we can easily apply such an argument to the 
two-point correlation of anisotropic fluctuations in the T 3 model. In that case, corresponds 
to an expansion coefficient a/ m and r is generalized to (ai m a*, m ,) / \jC\C\ for 1^1' orm^ m! . 

Thus, the presence of off-diagonal components in the covariance matrix causes a squeeze 
of the distribution function, leading to an increase in the best-fit normalization and a de- 
crease in the likelihood value. 
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